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Nonuniversal boundary conditions in grand unified theories can lead to nonuniversal gaugino 
masses at the unification scale. In i?-parity preserving theories the lightest supersymmetric particle 
is a natural candidate for the dark matter. The composition of the lightest neutralino and the 
identity of the next-to- lightest supersymmetric particle are studied, when nonuniversal gaugino 
masses come from representations of SO (10). In these cases, the thermal relic density compatible 
with the Wilkinson Microwave Anisotropy Probe observations is found. Relic densities are compared 
with the universal case. Mass spectra in the studied cases are discussed. 
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I. INTRODUCTION 

The phenomenology of supersymmetric models de- 
pends crucially on the compositions of neutralinos and 
charginos, if the lightest neutralino is the lightest super- 
symmetric particle (LSP). In addition to the laboratory 
studies, relevant input is obtained from the dark matter 
searches, where the the Wilkinson Microwave Anisotropy 
Probe (WMAP) satellite has put precise limits on the 
relic density. Supersymmetric theories which preserve 
i?-parity contain a natural candidate for the cold dark 
matter particle. Neutralino LSP can provide the appro- 
priate relic density. 

In many supergravity type models the lightest neu- 
tralino is binolike, which often leads to too large thermal 
relic density, as compared to the limits provided by the 
WMAP experiment [1]. When the gaugino masses are 
not universal at the grand unification scale, the resulting 
neutralino composition changes from the case of univer- 
sal gaugino masses [2]. In this paper, the thermal relic 
density of the neutralino LSP is studied, when gaugino 
masses are due to nonuniversal representations of SO (10) 
grand unified theory (GUT) [3, 4]. Dark matter in a 
particular gauge symmetry breaking chain of the SO (10) 
GUT in the case of universal gaugino masses has been 
recently studied in [5]. Some phenomenological aspects 
of SO (10) GUTs with nonuniversal gaugino masses have 
been considered in [6, 7]. 

SO (10) has many attractive features among the GUT 
models. One of the most appealing properties is that 
one family of matter fermions can be put into a sin- 
gle 16-dimensional irreducible spinor representation of 
SO (10), including the right-handed neutrino [8, 9]. In 



addition, SO (10) allows possibility for the Yukawa cou- 
pling unification and representations are anomaly free. 
Conservation of i?-parity, which forbids the unwanted 
dimension-five operators leading to rapid proton decay, 
may result from the SO (10) symmetry breaking. The 
doublet-triplet splitting could be achieved using, e.g. , 
the so-called Dimopoulos-Wilczek mechanism [10]. Be- 
cause the SO (10) gauge symmetry breaks down to the 
standard model (SM) gauge symmetry through some in- 
termediate group, the SO (10) GUT offers several possi- 
bilities for the model building. For example, it can con- 
tain as a subgroup the Pati-Salam SU(4) xSU(2) xSU(2) 
model. 

Gaugino masses originate from the non-renormalizable 
terms in the = 1 supergravity Lagrangian involving 
the gauge kinetic function fab{^) [H]- The gauge part of 
the Lagrangian contains the gauge kinetic function cou- 
pling with two field strength superfields W^. The La- 
grangian for the coupling can be written as 



(1) 



where a and b are gauge group indices (for example, 
a^b = 1,2, ...,45 for SO(IO)), and repeated indices are 
summed over. The function fab{^) is an analytic func- 
tion of the chiral superfields <l> in the theory. The chiral 
superfields <l> consist of a set of gauge singlet superfields 
<l>^ and gauge nonsinglet superfields <l>^ under the grand 
unified group. The gauge kinetic function fab{^) can be 
expanded, 



fabW 



(2) 
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where <l>^ and <l>^ are the singlet and nonsinglet chiral su- 
perfields, respectively. Here /o(*^*^) and fn{^^) are func- 
tions of gauge singlet superfields and Mp is some 
large scale. In order to generate a mass term for the gaug- 
inos, the gauge kinetic function must be non-minimal, 
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i.e., it must not be a constant [12]. When gets a 
vacuum expectation value {vev) {F^)^ the interaction (1) 
gives rise to gaugino masses: 

CgkD^-^^yx' + ^-C, (3) 

where A^'^ are gaugino fields. The nonuniversal gaugino 
masses are generated by the nonsinglet chiral superfield 
that appears linearly in the gauge kinetic function 
fabm in Eq. (2). 

Gauginos belong to the adjoint representation of the 
gauge group, which in the case of SO (10) is the 45 di- 
mensional representation. Because Eq. (3) must be gauge 
invariant, <l> and F^ must belong to some of the fol- 
lowing representations appearing in the symmetric prod- 
uct of the two 45 dimensional representations of SO (10) 
[13, 14]: 

(45 (g) 45)symm = 1 © 54 © 210 © 770. (4) 

The representations 54, 210 and 770 may lead to 
nonuniversal gaugino masses, while the 1 dimensional 
representation gives manifestly the universal gaugino 
masses. The relations between the gaugino masses are 
determined by the representation invariants, and are spe- 
cific for each of the representations. Because the gauge 
kinetic function in Eq. (2) can get contributions from 
several different <l>'s, a linear combination of any of the 
representations is also possible. In that case the gaug- 
ino mass terms are not uniquely determined anymore, in 
contrast to the contribution purely from one representa- 
tion. Here we assume that the dominant component of 
the gaugino masses comes from only one representation. 
This gives us a clear understanding of the role of different 
representations. 

II. DARK MATTER IN SO(IO) 
REPRESENTATIONS 

A. Breaking Chains: SO(IO) ^ H ^ SM 

The GUT group SO (10) breaks down to the standard 
model gauge group SU(3)xSU(2)xU(l) via some inter- 
mediate gauge group H. Therefore the gaugino mass 
relations depend also on the gauge group breaking chain, 
in addition to the representation invariants coming from 
the gauge kinetic function. Moreover, the intermediate 
breaking scale affects also the generated gaugino masses 
via heavy gauge supermultiplets that correspond to the 
broken generators. However, if the gauge breaking from 
the GUT group to the SM group takes place at the GUT 
scale, these loop-induced messenger contributions [15] 
can be neglected in comparison to the tree-level contri- 
butions. Some fits to the experimental data in SO (10) 
GUT indicate that the two breaking scales are very close 
to each other, see [16, 17], although realistic models ex- 
ist also with large splitting of the scales [18]. In this 



work we assume that the breaking from SO (10) to the 
SM gauge group happens at the GUT scale, and that the 
GUT breaking does not affect the gauge coupling unifi- 
cation. 

We will study the representations 54 and 210 in the 
right-hand side of Eq. (4). The interesting breaking 
chains of 54 and 210 are included also in the breaking 
chains of 770. Table I shows possible SO (10) breaking 
chains [14, 19, 20], which include the standard model 
gauge group, for the two chosen representations. Some 
of the subgroups lead to universal gaugino masses, or 
to massless gauginos [14], and we do not consider them. 
We will limit ourselves to the intermediate gauge groups 
SU(4)xSU(2)xSU(2), SU(2)xSO(7) and SU(5)xU(l). 



TABLE I: Breaking chains of SO (10) representations 
54 and 210 which include the SM gauge group. 



F^ 


H 


Subgroup description 




SU(4)xSU(2)xSU(2) 


Pati-Salam 


54 


SU(2)xSO(7) 






SO(9) 


Universal gauginos 


210 


SU(4)xSU(2)xSU(2) 
SU(3) xSU(2) xSU(2) xU(l) 
SU(3)xSU(2)xU(l)xU(l) 
SU(5)xU(l) 


Massless gluino 
Massless su(2)l gauginos 

"Flipped" su(5) 



Table II displays the ratios of resulting gaugino masses 
at the tree level as they arise when F^ belongs to the 
above-mentioned representations of SO (10) or singlet 
[14]. The resulting 1-loop relations at the electroweak 
scale are also displayed. These values and the resulting 



TABLE II: Ratios of the gaugino masses at the GUT 
scale in the normalization Ms{GUT) = 1, and at the 
electroweak scale in the normalization Ms{EW) = 1. 



F^ 


H 


Mf 












1 




1 


1 


1 


0.14 


0.29 


1 


54 


SU(4)xSU(2)xSU(2) 


-1 


-1.5 


1 


-0.15 


-0.44 


1 


54 


SU(2)xSO(7) 


1 


-7/3 


1 


0.15 


-0.68 


1 


210 


SU(5)xU(l) 


-96/25 


1 


1 


-0.56 


0.29 


1 



relic densities can be compared with the universal and 
nonuniversal representations resulting in the SU(5) GUT 
model [21, 22]. Since we assume breaking at one scale, 
the universal model with which we should compare in 
the SO (10) GUT is similar to the universal model in the 
SU(5) GUT. In the nonuniversal representations, the re- 
lations between gaugino masses change. Thus, e.g. , the 
54-dimensional Pati-Salam model of SO (10) may seem 
at first glance rather similar to 24 of SU(5), but we will 
see that the twice as large bino component has a large 
effect to the relic density. The bino and wino mass pa- 
rameters affect directly the lightest neutralino mass and 
properties. They also affect the value of the /i-parameter 
through the renormalization group equations (RGE) and 



3 



the radiative electroweak symmetry breaking (rEWSB), 
therefore controhing also the Higgsino component in the 
Ughtest neutrahno. Since the hghtest neutrahno mass 
hmit can be deduced from the chargino mass hmit, the 
nonuniversal gaugino masses change the lower limit for 
the neutrahno mass: for 54 the neutrahno mass limit is 
smaller than in the universal case, while for 210 the mass 
limit is close to the chargino mass limit. 

B. Calculation of Dark Matter Relic Density 

We calculate the SUSY spectrum for each model with 
the program SOFTSUSY (version 2.0.11) [23], and the 
resulting relic density with the program micrOMEGAs 
(version 2.0.7) [24, 25, 26]. For the relic density, we use 
here the WMAP combined 3 yr limits [1] 

Ocdm/i' = 0.1 1054™76 (2^). (5) 

In all the figures that we show below, the filling denoted 
by wmap is the WMAP-preferred region. For the b ^ 
experimental branching fraction, we have used the two 
sigma world average [27], 

BR{b 57) = (355 ± 24l?o ± 3) x 10"^. (6) 

The areas enclosed by the bsg contour are disallowed by 
the b ^ sj constraint. For the particle masses, the fol- 
lowing limits are applied [26]: mg^ > 99.4 or 100.5 GeV 
depending on if the lightest neutrahno mass is below or 
above 40 GeV, m^^ > 95 GeV, rrif^ > 80.5 to 88 GeV 
depending on the lightest neutrahno mass (from 10 to 75 
GeV), m^. > 43 GeV, and m^± > 73.1 to 103 GeV de- 
pending on the sneutrino masses (from 45 to 425 GeV). 
In the figures, lep shows an area where the experimental 
mass limits are not met, rge shows an area where there 
is no radiative EWSB, and Isp the area where neutrahno 
is not the LSP. The curve ruh = 114 GeV is depicted 
in the figures (dash-dotted line denoted by h). For the 
shown parameter regions, when otherwise experimentally 
allowed, Higgs is always heavier than 91 GeV, which is 
the Higgs mass limit in MSSM for tan/? > 10 assuming 
maximal top mixing [28]. 

1. Representation 54 

The area of preferred thermal relic density for the two 
chains of the 54 dimensional representation are shown for 
sets of parameters in Figs. 1, 2, and 3. In each set of three 
figures, the first figure (a) represents the neutrahno relic 
density for given parameters with collider constraints de- 
picted in the plot, the second figure (b) shows, for the 
same parameters, the identity of the next-to-lightest su- 
persymmetric particle, and the third figure (c) shows the 
lightest neutrahno composition in RGB-color encoding, 
{i.e., colors, or hues of black and white, indicate the par- 
ticle as shown in the figure; therefore the mixture of the 



colors, or hues of black and white, describes the nature of 
the Xi-composition). In each of the figures, the WMAP- 
preferred relic density filling is also superimposed to the 
graph. 

As can be seen from the Table II, the lightest neu- 
trahno is expected to be bino rather than wino. The 
large bino component tends to suppress the neutrahno 
annihilation cross section, since bino lacks the s-channel 
Z-boson annihilation mode. A substantial Higgsino com- 
ponent is usually needed to help to increase the annihi- 
lation rate, unless there happens to be coannihilation or, 
e.g. , an open Higgs s-channel annihilation mode avail- 
able. 

In Fig. 1 relic density, the next-to-lightest supersym- 
metric particle (NLSP) and LSP composition in the 
breaking chain SU(2)xSO(7) are shown. Because Xi is 
mostly bino, the spectrum with preferred relic density is 
quite light and confiicts with collider constraints in some 
parts of the parameter space. With increasing gaugino 
masses also the Higgsino component in the neutrahno 
LSP increases, and at the point where the change to dom- 
inantly Higgsino LSP occurs, also the relic density drops. 
The overall relic density is not very high, thus allowing a 
wider WMAP-preferred region than, e.g. , in the singlet, 
i.e., mSUGRA case [22]. For a given M2, the correspond- 
ing Ml is smaller than in the singlet case, which results 
in a smaller /a value at the EW-scale. This has an effect 
of an increasing the Higgsino component in the lightest 
neutrahno thus boosting the annihilation. The NLSP 
is chargino, and with increasing Higgsino component it 
eventually becomes the LSP. 

In the allowed region, M2 is less than 740 GeV, which 
restricts the lightest chargino mass to values less than ~ 
150 GeV. The lower limit for the chargino mass is the 
LEP limit. The partners of the SM fermions are heavier 
than 300-500 GeV. Thus, assuming that the neutrahno 
is responsible for the dark matter, this breaking chain 
of the gauge symmetry has as a robust prediction for 
the upper limit of the chargino mass, and furthermore 
it is lighter than the squarks and sleptons. For M2 ~ 
350 GeV the whole neutrahno and chargino spectrum, 
and even the gluino, is lighter than the sfermions for the 
LEP allowed and WMAP-preferred region. The same is 
true also for M2 - 630 GeV and mo > 940 GeV for the 
WMAP-preferred region. 

The b ^ sj constraint cuts away a considerable area 
from the otherwise allowed region. Including a 10% error 
in theoretical calculations of the decay, the constraints 
from 6 ^ 57 loosen considerably, and all of the other- 
wise allowed WMAP area becomes available. The lightest 
Higgs is always lighter than 114 GeV, but heavier than 
91 GeV. 

The effect of varying the universal trilinear coupling Aq 
is shown in Fig. 2. In contrast to the previous figure, here 
the sign of the /i-parameter is chosen to be negative. In 
general, giving the Ao-parameter a nonzero value tends 
to increase the relic density. Increasing \Aq\ will help 
the m'jj to run larger negative values during the KG- 
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FIG. 1: Relic density Q^h'^ in the representation 54 with H =SU(2) xS0(7) in the (M2,mo) plane for 
tan/3 = 10, sgn(/i) = +1, Aq = 0. In la the dark shaded areas represent the larger relic density. The filling denoted 
by wmap is the WMAP preferred region, lep shows an area, where the experimental mass limits are not met, rge 
shows an area where there is no radiative EWSB, and Isp the area where neutralino is not the LSP. h encloses the 
area with rrih < 114 GeV, and in the following figures bsg the area disallowed hy b ^ limits. In lb the NLSP 
within the same region is plotted, and in Ic the Xi composition. 



evolution, and therefore to increase the actual (absolute) 
value for the /i-parameter via the rEWSB. This in turn 
favors the bino component in the Xi composition in the 
SU(2)xSO(7) chain. The effect can be seen in the LSP- 
composition Fig. 2c. The preferred WMAP region still 
follows the transition zone of the Xi from bino to Hig- 
gsino. The effect of negative /i is most visible in the fact 
that the NSLP in this case is X27 since negative /i tends to 
increase the lightest chargino mass (at least in the limit 
of IM2I < The b ^ constrains an area with neg- 
ative Aq values at M2 < 1 TeV, and again including also 
the theoretical error, the constrained area becomes con- 
siderably smaller. The light Higgs boson mass is heavier 
than 114 GeV above the dash-dotted line. 

In Fig. 3 the relic density, NLSP and x^ composition 
are plotted for the SU(4) xSU(2) xSU(2) breaking chain. 
In this breaking chain the | Mi | and | M2 \ values are closer 
to each other at the EW-scale than in the SU(2)xSO(7) 
chain, but still wider spread than in the singlet. That 
has an effect of increasing the /i value, therefore result- 
ing in a smaller Higgsino component for the equal M2 
values for the two chains of the 54 dimensional represen- 
tation. Again the WMAP allowed narrow region follows 
the transition from the Xi from bino to Higgsino. The 
b ^ and Higgs boson 114 GeV limits cut pieces from 
near the low M2, mo-values. The NLSP is mostly the 
lightest chargino, but an interesting region exists with 
small mo near the area, where Xi is no longer the LSP; 
along the line of transition from the stau NLSP to the 
smuon NLSP there is a narrow region, where Xi, stau and 
smuon masses are very close to each other, and the coan- 
nihilations may reduce the relic density to an acceptable 
level. For example, for M2 = 1400 GeV, the LSP mass is 



around 410 GeV, and for M2 = 1000 GeV, the LSP mass 
is 290 GeV for the WMAP-preferred area. However, this 
area is highly prone to numerical subtleties, and the or- 
dering of the LSP identity changes in the preferred relic 
density range, when comparing the output of different 
spectrum calculators. Reducing the mo parameter fur- 
ther makes the stau become the LSP. From the collider 
point of view, such regions may be especially interesting, 
as they would lead to quasistable smuons or staus. 



2. Representation 210 

In the representation 210 we inspected the break- 
ing chain through the intermediate gauge group 
SU(5)xU(l), called fiipped SU(5) [29, 30, 31, 32]. In 
Fig. 4a the area of preferred thermal relic density in the 
representation 210 is plotted for a set of (GUT scale) 
parameters. For the chosen parameters, rather large 
WMAP-preferred regions are found for large values of M2 
and/or mo parameters. When not Higgsino, the lightest 
neutralino is expected to be wino, rather than bino (see 
Table II), and therefore the neutralino relic density to 
be very small. In general, due to the wino being the 
smallest of the two electroweak gaugino parameters, it 
characterizes the lightest neutralino. Since the lightest 
chargino is characterized also by this parameter, for a 
large part of the parameter space, the masses of the Xi 
and xt very close to each other, which boosts the 
rapid neutralino relic density annihilation. The situation 
resembles the one arising in the anomaly mediated SUSY 
breaking scenario, where also both the lightest neutralino 
and chargino are characterized by the wino mass param- 
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FIG. 2: Representation 54: H =SU(2) xS0(7), tan^^ = 10, sgn(/i) 
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FIG. 3: Representation 54 with H =SU(4) xSU(2) xSU(2) in the (M2,mo) plane for 
tan^^ = 10, sgn(/i) = +1, Aq = 0. Otherwise as in Fig. 1. 



eter. 

Since both the wino and Higgsino have a large anni- 
hilation cross section, the WMAP-preferred relic density 
region does not have to follow the transition zone of Xi 
from wino to Higgsino. The increase of the relic den- 
sity to the observed level is mainly due to the increase of 
the mass parameters M2 and mo. An interesting change 
in the pattern can be seen on the diagonal of the fig- 
ures, where the Xi and Xi •> ^"^^ also at some point X2-> 
masses are very close to each other leading to enhanced 
coannihilation through the processes Xixf ^ QuQd and 
XiXi^xfxf qq^i£^W~^W~ ^ which allows the accept- 
able parameter space to extend to very large values of 
M2 and mo. 

The WMAP-preferred region is very wide as compared 
to, e.g. , the universal model or other models in SU(5) 
[22, 33]. The spectrum is relatively heavy for the WMAP- 
preferred region, around a couple of TeV. This leads to 



a wide range in the parameter space, but from the point 
of view of the coming Large Hadron Collider such a mass 
spectrum may be problematic. However, the represen- 
tation 210 produces naturally a neutralino with mass 
around a TeV, which seems to be favorable in the view 
of the recent PAMELA [34] and ATIC [35] results of the 
excess positron and positron+electron flux. Under cer- 
tain circumstances, a nearby clump of 600-1000 GeV neu- 
tralino LSP could fit into these observations [36]. 

In Fig. 5 the area of preferred thermal relic density 
in the representation 210 is plotted for the same set of 
(GUT scale) parameters as in Fig. 4, except that now 
the trilinear Ao-parameter is varied along the y-axis, and 
the mo is set to 1 TeV. The interesting feature in this 
figure is the existence of the pseudoscalar Higgs anni- 
hilation channel through the M2-values. This reduces 
greatly the relic density in the parameter space where the 
LSP mass equals or is less than half of the A-Higgs mass. 
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FIG. 4: Representation 210: tan^^ = 10, sgn(/i) = +1, Aq = 0. Otherwise as in Fig. 1. 
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This has the effect of pushing the WMAP-preferred rehc 
density region to heavier neutrahno masses, and there- 
fore to larger M2 values. The width of the WMAP re- 
gion is naturally the same as before, since the top of 
the figure, where = 0, coincides with Fig. 4 with 
mo = 1000 GeV. With these parameters, the NLSP is 
always the lighter chargino, and the lightest neutralino is 
dominantly a wino, which can also be read from Fig. 4c. 



In the representation 54, the lightest neutralino is 
predominantly a bino, leading to the narrow areas of 
the WMAP-favored region. The excessive relic density 
is diluted either by increasing the Higgsino component 
or by coannihilation with other particles. The break- 
ing chain SU(2)xSO(7) predicts an upper limit for the 
lighter chargino mass for the chosen parameters. For 
part of the WMAP allowed region, the whole neutralino 
and chargino spectrum is lighter than the spectrum of 
sfermions. For the SU(4)xSU(2)xSU(2) breaking chain, 
the relic density area is narrow in the parameter space. 
Interestingly, there may exist a region, where the stau, 
smuon and the lightest neutralino masses are in a very 
close range to each other. This can lead to long-lived 
staus and smuons, which may be stable in the collider 
time scale. 



In the 210 dimensional representation the lightest neu- 
tralino is either wino or Higgsino, which leads to a low 
thermal relic density. In addition, the lightest chargino 
and the lightest neutralino tend to be close in mass, 
thus providing a coannihilation channel. The preferred 
relic density area is quite large. The sparticle spectrum 
is heavy, as compared to the universal mSUGRA case. 
Only in a small part of the WMAP-preferred parameter 
region are a few SUSY particles expected to be within the 
kinematic reach of the LHC. 



III. DISCUSSION AND SUMMARY 

We studied the dark matter allowed regions in the 
SO (10) GUT representations, of which all but the singlet 
may lead to nonuniversal gaugino masses. The WMAP- 
preferred relic density regions are quite distinct for differ- 
ent representations, thus leading to quite different parti- 
cle spectra for each representation. 



In this work we have studied each representation sepa- 
rately. It is obvious that if several representations affect 
simultaneously the composition of neutralinos, the possi- 
ble WMAP-preferred region in the parameter space may 
be relaxed. However, if the neutralinos and charginos are 
found with a certain mass pattern, it helps to understand 
the relation of the lightest neutralino with dark matter, 
if characteristics of each representation are known. 
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